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Abstract 

A brief  outline  of  the  mathematical  technique  of  scalar  Hertz  potentials  is  provided.  It  is 
shown  that  the  technique  can  be  successfully  implemented  if  all  four  constitutive  dyadics  of 
a linear  bianisotropic  medium  have  a gyrotropic  structure.  A medium  of  that  nature  is  a 
so-called  Faraday  chiral  medium,  and  it  may  even  be  nonhomogeneous  but  of  a restricted 
form. 

1.  Introduction 

Approaches  to  solve  electromagnetic  field  problems  can  be  roughly  divided  into  two  categories: 
direct  and  indirect.  While  the  former  deal  directly  with  the  Maxwell  equations  for  the  electro- 
magnetic field  vectors,  the  latter  aim  to  derive,  and  consequently  solve,  alternative  differential 
equations  that  are  fulfilled  by  certain  functions  — scalar,  vector  or  dyadic.  The  field  vectors 
themselves  then  follow  from  these  functions  through  (relatively  simple)  manipulations  such  as 
differentiations  and  integrations. 

Two  of  the  most  prominent  indirect  solution  techniques  are  (i)  Green  function  methods 
and  (ii)  potential  methods.  In  this  contribution  we  shall  focus  on  potentials,  and  specifically 
on  scalar  Hertz  potentials  and  their  role  in  the  solution  of  electromagnetic  field  problems  in 
complex  mediums. 

Vector,  scalar  and  Hertz  potentials,  in  their  application  to  (achiral)  isotropic  mediums,  have 
been  prominent  mathematical  tools  since  the  early  stages  of  the  development  of  the  Maxwellian 
formalism  — a fact  easily  betrayed  by  their  inclusion  in  the  vast  majority  of  standard  textbooks 
on  electromagnetic  theory.  But  only  in  the  last  two  decades  has  it  been  possible  to  extend 
the  concept  of  scalar  Hertz  potentials  into  the  realm  of  more  complex  mediums,  which  may  be 
anisotropic,  bianisotropic,  nonhomogeneous,  in  their  nature. 

Here,  an  outline  of  the  scalar  Hertz  potential  method  is  given,  necessarily  short  due  to  the 
limited  available  space.  Full  details  including  a comprehensive  listing  of  references  of  the  relevant 
research  literature  can  be  found  in  a recent  book  chapter  [1]. 

2.  Analysis 

The  Maxwell  equations  for  a general,  linear  bianisotropic  medium  are1 

ioj  [g  • E(x)  + £ • H(x)j  + Vx  H(x)  = Je(x) , (1) 

V x E(x)  - iw  [C  * E(x)  + fx  • H(x)j  = -Jm(x) . (2) 

‘Vectors  are  in  bold  face  whereas  dyadics  axe  in  normal  face  and  underlined  twice.  Contraction  of  indices  is 
symbolized  by  a dot  and  ab  is  a dyadic  product;  the  unit  dyadic  is  /. 
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E(x)  and  H(x)  are  the  electric  and  magnetic  field  vectors,  whereas  Je(x)  and  Jm(x)  are  the 
prescribed  electric  and  magnetic  current  densities  (a  time-dependence  of  exp(— iujt)  is  implicit). 
At  this  stage  no  restrictions  are  imposed  on  the  specific  form  of  the  four  constitutive  dyadics  in 
(1),  (2). 

To  develop  a field  representation  in  terms  of  scalar  Hertz  potentials,  a partial  scalarization 
of  the  differential  equations  (1),  (2),  is  carried  out  with  respect  to  an  arbitrarily  chosen  direction 
specified  by  a unit  vector  c — once  again,  full  mathematical  details  can  be  found  in  [1].  The 
decompositions 

E = Et  + Ecc,  H = Ht  + tfcc,  (3) 


are  introduced,  whereby  Et  • c = 0 and  Ht  • c = 0;  E<  and  Hf  being  called  the  transverse , Ec 
and  Hc  the  longitudinal  components  of  the  field  vectors.  Similarly,  for  the  electric  and  magnetic 
current  density: 

J e = Jet  "b  Jec  J m — J mt  "b  Jmc  ® > (4) 

and  also  for  the  derivative  operations: 


V = V, 


d 

+ C dxc  ' 


V2  = V?  + 


dx\ 


(5) 


wherein  V2  is  the  Laplace  operator. 

The  c-components  of  (1),  (2)  provide  two  scalar  expression  for  the  components  Ec  and  Hc 
and  one  finds,  after  some  algebra, 


Ec  = P ee  * Et  + P eh  • Ht  + Ke  , (6) 

Hc  = P he  * E*  + P hh  • H<  + Kh  , (7) 

where  Pee,  P eh,  P he  and  P hh  are  2x2  dyadic  operators  while  Ke  and  Kh  are  source  terms. 
Now,  a system  for  the  transverse  components  Et  and  Ht  is  derived  in  the  form 

3E 

_&T'b*a;=ee  *Et  + *wMeh  *Ht  = qte,  (8) 

3H 

~!hrc  + iUJ=he  * E‘  + i0JM=hh  ' H‘  = • (9) 

where  M.ee , M.eh,  M-he  and  M.hh  are  2x2  dyadic  differential  operators  of  second  order  and  qte 
and  qt/i  are  2- vector  source  terms. 

The  crucial  step  is  now  the  introduction  of  scalar  potentials  through 

Et  = Vt  $ + Vt  x 0c,  (10) 

Ht  = Vtn  + Vt  X ¥c.  (11) 

Straightforward  manipulation  of  (6),  (7),  then  leads  to  formulas  for  the  longitudinal  components 
Ec  and  Hc,  in  symbolic  form: 


Ec  — Ec  (<& , © , II , ^ ; Jec , Jmc)  j (12) 

tfc  = i?c($,0,n,$;Jec,Jmc)  . (13) 

The  expressions  (10)— (13)  constitute  a representation  of  all  field  components  in  terms  of  the  four 
scalar  functions  <&,  ©,  fl  and  'If;  substitution  into  the  available  differential  equations  produces  a 
system  of  four  differential  equations  for  the  four  scalar  potentials.  At  that  stage  no  direct  gain 
has  been  achieved  as  one  has  simply  exchanged  solving  a system  of  partial  differential  equations 
of  second  order  for  Et  and  Hf  with  solving  a corresponding  system  for  the  four  scalar  potentials. 
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The  question  is  therefore:  can  two  of  the  four  scalar  potentials  be  eliminated  in  a simple 
manner  so  that  the  field  vectors  E and  H can  be  represented  in  full  by  two  scalar  Hertz  potentials 
plus,  perhaps,  some  auxiliary  functions  that  are  required  to  deal  with  specific  source  terms? 

The  answer  to  that  question  is  that  such  an  elimination  is  not  possible  for  the  general, 
linear  bianisotropic  medium  characterized  by  e,  £,  £ and  /it.  Conditions  on  the  structure  of 
these  constitutive  dyadics  must  be  imposed  and  The  desired  reduction  is  possible  only  if  the 
constitutive  dyadics  are  of  the  form 


a = 


otL+  (o"c  — ct)  c c — iog  cx  l 


( °t 


~ia9 

V 0 


i<?n  0 \ 

<*t  o , 

0 (Jc  J 


£ = |,£,C,M- 


(14) 


All  four  constitutive  dyadics  are  of  gyrotropic  form,  and  a medium  of  that  nature  is  called 
a Faraday  chiral  medium  [2].  In  the  dyadics  (14),  at  and  ac  characterize  the  uniaxiality  of 
the  medium  whereas  ag  is  responsible  for  the  gyrotropic  or  gyrotropic-like  nature.  The  scalar 
Hertz  potential  formalism  for  homogeneous  Faraday  chiral  mediums  has  been  outlined  in  detail 
elsewhere  [2].  The  field  representation,  for  all  the  types  of  mediums  for  which  the  scalar  Hertz 
potential  technique  is  then  applicable,  can  symbolically  be  given  by 


E E (0  , ’F  J Jec  , Ue  , Ve  , Jmc  i um  > vm ) > (15) 

H H (0  , 'F  ; Jec  ) Ue  > ve  i Jmc  i um  i vm)  > (16) 


because  usually  the  two  potentials  that  get  eliminated  are  $ and  n. 

Typically,  the  two  functions  0 and  'F  — the  scalar  Hertz  potentials  — fulfil  a system  of 
differential  equations  of  the  structure 


Li©  + L2V  = si , (17) 

L3@  + L^  = s2,  (18) 

where  the  exact  form  of  the  second  order  partial  differential  operators  Li,  L2,  X3  and  L4  depends 
on  the  specific  type  of  bianisotropic  medium  being  analysed;  si  and  s2  are  source  terms.  In 
general,  thus,  0 and  'F  are  coupled  but  decoupling  occurs  if  the  medium  becomes  sufficiently 
simple  (as  happens,  for  example,  for  a uniaxial  dielectric  medium  [3]). 

Finally,  what  is  the  significance  of  the  functions  ue,  um,  ve  and  vm  that  have  entered  the 
field  representation  (15),  (16)?  These  are  auxiliary  functions  defined  by 

Jet  = Vtue  + Vt  X Vec,  (19) 

Jmt  = Vt  Um  + V(  X vm  c . (20) 

Consequently,  the  four  auxiliary  functions  are  calculated  from  the  differential  equations 

V^up  = Vt*Jpt,  V^Up  = Vt  • (c  x Jpt)  , p = e,m,  (21) 

and  it  is  immediately  clear  from  these  expressions  that  they  have  nothing  to  do  with  the  medium 
as  such  but  are  simply  generated  by  the  mathematical  formalism  to  accommodate  transverse 
current  densities. 

3.  Discussion 

The  brief  outline  of  the  scalar  Hertz  potential  technique  is  concluded  here  with  a number  of 
observations;  see  [1]  for  more  detailed  discussions  of  the  individual  points. 
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• Not  mentioned  previously  is  the  possible  generalization  to  nonhomogeneous  mediums, 
where  the  nonhomogeneity  is  of  a restricted  kind.  Namely,  all  constitutive  parameters 
can  depend  on  xc(=  x • c),  the  coordinate  parallel  to  the  distinguished  axis  c.  Thus,  the 
most  general  type  of  medium  that  is  amenable  to  the  scalar  Hertz  potential  technique  is 
characterized  by  constitutive  dyadics  of  the  form 


g(xc)  = et(xc)l  + 

ec(xc)  - et(xc)  ] 

cc  — ieg(xc)  c x I , 

(22) 

|(zc)  =(t(xc)L+  | 

£c(zc)  - Zt(xc) 

cc  - i£g(xc)  c x I, 

(23) 

C{xc)  = Ct(xc)L  + 

Cc(zc)  - (t{xc) 

] c c - i(g(xc)c  X X, 

(24) 

fjfxc)  = Pt(xc)L  + 

Pc(xc)  ~ Pt(xc)  c c — ipg(xc)  c x £. 

(25) 

All  constitutive  dyadics  have  a gyrotropic  structure  and  are  spatially  nonhomogeneous 
through  the  dependence  on  xc  and  are  thus  nonhomogeneous  Faraday  chiral  mediums  [2]. 

• The  auxiliary  functions  ue,  um,  ve  and  vm  are  determined  by  (19)— (21);  it  is  noted  that 

Jet  = 0 =4>  «e  = 0,  We=0,  (26) 

Jmt  = 0 k um  — 0 > vm  — 0 • (27) 

Consequently,  they  can  be  completely  omitted  from  the  mathematical  formalism  if  the 
medium  is  free  of  sources,  i.e.,  Je  = 0,  Jm  = 0,  or  if  the  current  densities  are  purely 
longitudinal,  i.e.,  J e = Jec c,  Jm  = Jmc c.  « 

• With  the  fundamental  system  of  differential  equations  (that  for  the  scalar  Hertz  potentials) 
that  needs  to  be  solved  being  a scalar  problem,  Green  functions  that  may  be  introduced 
subsequently  can  be  of  the  simpler  scalar  rather  than  much  more  complicated  dyadic 
nature. 

• The  introduction  of  the  scalar  Hertz  potentials  is  not  a unique  process.  New  scalar  Hertz 
potentials  0neu,(x)  and  'hneu;(x)  may  be  defined  via 

©netu(x)  = 7x  0(x)  + 72  'F(x) , ^netu(x)  = 73©(x)  +74  'F(x) , (28) 

where  7„,  n = 1,2, 3,4,  are  arbitrary  complex  constants  (provided  7174  7^  7273). 

In  conclusion,  the  scalar  Hertz  potential  technique  provides  a very  successful  mathematical 
method  to  deal  with  electromagnetic  field  problems  in  complex  mediums. 
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